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Abstract. Fayans energy density functional (EDF) FaNDF 0 has been applied to the nuclei around uranium 
region. Ground state characteristics of the Th, U and Pu isotopic chains, up to the two-neutron drip line, 
are found and compared with predictions from several Skyrme EDFs. The two-neutron drip line is found for 
FaNDF 0 , SLy4 and SkM* EDFs for a set of elements with even proton number, from Pb up to Fm. 


1 Introduction 

Presently, the Hartree-Fock (HF) or Hartree-Fock-Bogo- 
liubov (HFB) methods together with the effective Skyrme 
forces [|I1], Gogny forces or relativistic mean-field 
(RMF) models yj] are most common microscopical mod¬ 
els applied for description of the ground state properties 
of the heavy nuclei. All of the aforementioned approaches 
are usually interpreted as a version of the energy density 
functional (EDF) method suggested by Kohn and Sham 
Hi . This method is based on the theorem of Hohen- 
berg - Kohn yfl, which states that the ground state en¬ 
ergy Eq of any quantum system is a functional of its den¬ 
sity p( r). By itself, the theorem says nothing about the 
form of this functional, and various options for the Skyrme 
EDF and the Gogny EDF are, in fact, different “ansatzs”. 
Among popular Skyrme EDFs, there are quite old func¬ 
tionals SKM* @] and SLy4 0, see the review article 0. 
A record in the accuracy of the description of the nuclear 
masses belongs to the Skyrme functional HFB-17 [J9|] with 
the average deviation from the experiment being around 
600 keV, obtained by adding phenomenological correc¬ 
tions atop of the mean-field. Of a comparable accuracy are 
the other functionals of this family, up to the HFB-27, the 
corresponding tables are presented on the site 11Q0- Also 
relatively new functional UNEDF1 0 has proved to be 
very successful in describing the deformed nuclei. 

It is also worth to mention a relatively newly developed 
approach, known originally as a BCP (Barcelona - Cata¬ 
nia - Paris) [12] method, and later as a BCPM (Barcelona 
- Catania - Paris - Madrid) method 0. The main bulk 
term of the BCPM functional was found by starting from 
the equation of the state of nuclear and neutron matter. 


obtained within the Brueckner-Hartree-Fock method by 
using a realistic NN + NNN potential. 

We use the EDF developed by S.A. Fayans with coau¬ 
thors [14-lTt]. In comparison to the Skyrme or Gogny 
EDFs, it possess two main peculiarities. Firstly, the main 
in-volume term of the Fayans EDF has more sophisticated 
form. It can be schematically written as 


E(p) = 


ap 2 1 + ap° 
2 1 + yp a 


( 1 ) 


a e-mail: tolkn@mail.ru 
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The corresponding term of the Skyrme EDF would cor¬ 
respond to y -0 case within this relation. The use of the 
bare mass, m*=m, is another peculiarity of the Fayans 
method. Both features of this approach are closely re¬ 
lated to the self-consistent Theory of Finite Fermi Sys¬ 
tems (TFFS) [181. The latter is based on general princi¬ 
ples of the TFFS 1 19 1 sup plemented with the TFFS self- 
consistency relations [20]. These two peculiarities of the 
Fayans approach reflect, in a effective manner, energy de¬ 
pendence effects inherent to the self-consistent TFFS. E.g., 
the effective mass in this approach is a product m*=m* k ■ m* E 
of the “k-mass” and the “E-mass”. The two effects com¬ 
pensate each other almost exactly 0 resulting in m* - 1. 
In the Skyrme HF method, the E-mass is identically equal 
to unit and the effective mass may deviate from unity. 

Until recently, the Fayans method was applied for 
spherical nuclei only. It turned out to be rather success¬ 
ful in systematic description of nuclear magneticJ21. 22] 
and quadrupole [23l!24|l moments, nuclear radii l25l l. beta- 
decay probabilities [26], the energies and B(E2) values for 
the first excited 2 + states in semi-magic nuclei [27| 28]. 
In a recent study of the single-particle energies in seven 
magic nuclei, a record accuracy was achieved [2911 . Most 
of these calculations were made with the EDF DF3-a 0. 
which is a small modification of the DF3 00 param- 
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eter set, concerning the spin-dependent terms of the EDF. 
In particular, the effective tensor term of DF3-a is signifi¬ 
cantly stronger than the DF3 one. 

The Fayans EDF was applied recently to deformed nu¬ 
clei for the first time 0. A localized version FaNDF 0 
0 of the general finite range Fayans EDF was used 
which makes its surface term more similar to the Skyrme 
one. This allowed to employ the computer code HF- 
BTHO l32ll . developed for Skyrme EDFs, with some mod¬ 


ifications. 

The results obtained with the FaNDF 0 EDF for Pb 
and U isotopic chains turned out to be promising. In 
the present work, we continue these calculations around 
the uranium region. In particular, we investigate the two- 
neutron drip line for a set of elements from Pb up to Fm, 
those with even proton number being considered. For a 
systematic comparison, we carried out also all calculations 
for two popular Skyrme EDFs, SkM* [@] and SLy4 0]. 
For completeness, we included also into the analysis pre¬ 
dictions from the HFB-17 and HFB-27 functionals taken 
from Ref. [1(J. In addition we discuss deformation char¬ 
acteristics of the drip line nuclei. 

Recent interest to the problem of fixing the neutron 
drip line [ 33. i34] is partially induced with importance of 
this characteristic of the nuclear chart for analysis of the 
r-process dynamics in stars |35j-37]. A couple of remarks 
should be made concerning validity of the EDF method 
with fixed set of parameters for predicting the drip lines. In 
particularly for the older generation EDFs, the parameter- 
izations were obtained by using data on nuclei close to the 
stability, that is, with the chemical potential-8 MeV. 
When extrapolating to very neutron rich systems, the role 
of EDF parameter errors becomes prominent 1331 38]. 
Also, analysis of the problem on the basis of the many- 
body theory point of view |[39| |40] shows that in vicinity 
of drip lines, that is p n —» 0 or p p —* 0, the EDF parame¬ 
ters describing the effective A'A'-interaction at the nuclear 
surface may need to be modified. Close to the neutron 
drip line, with a small chemical potential p n , attractive part 
of the /V/V-interaction may become enhanced, resulting a 
deeper neutron mean field. This, as a consequence, shifts 
the drip line farther away. To take this effect into account, a 



Figure 1. S^n values in the uranium chain for different EDFs. 
Experimental and estimated values are taken from ITfl l. 



Figure 2. /J 2 values in the uranium chain for different EDFs. 


simple model was used in [39|,|40] which, evidently, over¬ 


estimates the effect, predicting unrealistically strong shift 
of the drip line. It was, nevertheless, demonstrated that 
this effect shifts the drip line. However, one should bear 
in mind that by explicitly introducing such kind of compo¬ 
nent to the EDF model, the model parameters need to be 
readjusted, which, so far, has not yet been done. Thus, the 
total net effect remains still unknown. In the present work 
we use the standard approach with fixed EDF parameters. 


2 Deformation properties of Th, U, and Pu 
isotopic chains 

Details of the computation scheme are identical to those 
described in 113ill . We employ the axial computer code 
with the oscillator basis, the number of the oscillator 
shells being equal to N s h=25. We limit ourselves to the 
quadrupole deformation /?2 only, with reflection symmetry 
assumed. All the parameters of the normal component of 
the used FaNDF 0 EDF are the same as in 0. As for the 
anomalous term of the EDF 0, 

£anom = £ vJ(r)F^ (p + (r))v,(r), (2) 


where v,(r) is the anomalous density, a simplified version 
is used, 

F?(p+) = C 0 (£ + h ( (p+/po)) . (3) 

Here p+-p„+p p , and C'o-n 2 /ppm is the usual for TFFS 
normalization factor. The HFB equations are solved by 
using quasi-particle cut-off energy E cul = 60 MeV. The 
main part of calculations are carried out by employing 
simplest (“volume”) model of pairing, If = 0, with 
f£ x = -0.440. All calculations are repeated for two 
Skyrme EDFs, SKM* and SLy4, the results being com¬ 
pared also with predictions [ 10] from another two Skyrme 
EDFs, HFB-17 and HFB-27 EDFs. 

Let us begin from the uranium chain. The two-neutron 
separation energies. 


S 2n (N, Z ) = B(N, Z) - B(N - 2, Z), (4) 

are displayed in Fig. Q] Here B(N,Z) is the binding en¬ 
ergy of the nucleus under consideration. Comparison is 
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Figure 3. £ def values in the uranium chain for different EDFs. 


Figure 4. S 2 „ values in the thorium chain for different EDFs. 


made with experimental data [41] and predictions from 
four Skyrme EDFs. We first consider region of A < 244 
with known experimental values. The HFB-17 and HFB- 
27 models reaffirm their high accuracy. As to the FaNDF 0 
functional, agreement also looks rather reasonable, taking 
into account that the parameters were fitted d only for 
spherical nuclei not heavier than lead. The deviation of 
0.5 MeV from the experimental S 2 ,, values for heavy U 
isotopes is explained mainly by two reasons, with is the 
use of a simple volume pairing interaction, and absence 
of the effective tensor term in the FaNDF 0 EDF. Indeed, 
as it was shown in |[30], the tensor term is especially im¬ 
portant in uranium and transuranium region as, in corre¬ 
sponding spherical case, high -j levels dominate in vicin¬ 
ity of the Fermi level for these nuclei. As a result, the 
spin-orbit density, which comes to the EDF together with 
the tensor force, is typically large in these nuclei, chang¬ 
ing significantly along the isotopic chain. For the SFy4 
EDF agreement is a bit worse, whereas the disagreement is 
more significant in the SKM* case. It is worth to note that 
this Skyrme EDF, being fitted to masses not with so high 
accuracy as some modern equivalents, reproduces e.g. fis¬ 
sion properties in actinides or energies and B(E2) values 
of the first excited 2 + states in semi-magic nuclei [42] rel¬ 
atively well. For higher A values the SKM* values of S i„ 
are significantly higher than those for all other EDFs. As 
a result the corresponding drip line point 4 < 1 l '' p -320 turns 
out be significantly farther away than for all other EDFs. 
This quantity is defined as the last nucleus for which the 
two-neutron separation energy is yet positive. For SFy4 
we note that S 2 n at A - 286 is slightly negative. However, 
after this point the S 2 n remain positive up to A = 298. This 
kind binding re-entrance was also predicted at l(33ll . Thus, 
we can set A^CSFyT) = 298. 

The ground state quadrupole deformation parameter /T 
of the U chain and the corresponding deformation energy. 


£dcr(yS 2 ) = B(fr) - B(p 2 = 0), (5) 

are displayed in Figs. [2]and[3] correspondingly. Here, each 
curve is cut around the corresponding drip point. Gener¬ 
ally, all five EDFs under discussion behave in similar way, 
the most of nuclei in this chain being deformed in a pro¬ 
late way whereas there exists a region of spherical nuclei 


around ,4-280. In more detail, the width of the spherical 
region is narrowest for the FaNDF 0 EDF and widest for 
SFy4. Both of the HFB EDFs predict oblate deformations 
for four nuclei 262<A<268, whereas the FaNDF 0 EDF, 
only for two of them. A -264 and A— 266. 

The reason for this difference is quite simple. In vicin¬ 
ity of the phase transition, with a change of the deforma¬ 
tion sign, there are typically two energetically close by 
minima, the prolate and oblate one. Their closeness may 
be confirmed with the observation that the deformation en¬ 
ergy, see Fig. |3] does not show any non-regular behav¬ 
ior in the transition point, where the order of two minima 
changes. In such kind of situation, the transition value of 
A may move even due to a small change of the EDF pa¬ 
rameters. In the mass region of A > 280 the prolate defor¬ 
mation arises for the FaNDF 0 EDF for three points earlier 
than for both HFB and SFy4 EDFs. In this region the max¬ 
imum value of /T-0.25 appears for FaNDF 0 EDF just in 
the drip line point. The corresponding value is a bit less 
for the HFB-17 EDF, being /T-0.2 for other three EDFs. 
The difference is greater for the deformation energy due 
to its quadratic behavior around the minima. Thus, all five 
EDFs under consideration predict a well developed prolate 
deformation for uranium isotopes in the vicinity of the drip 
line. 

Next, we investigate the thorium chain. The two- 
neutron separation energies are displayed in Fig. |4] Again, 
all EDFs, with exception of SKM*, behave in general sim¬ 
ilarly in the region of A < 280, the latter curve being sig¬ 
nificantly higher. Close to the drip line, however, the pic¬ 
ture is different. In accordance with the above discussion, 
binding re-entrance can be now seen also for FaNDF 0 . As 
a result, we obtain A" p (FaNDF°)=294. For the HFB-17 
and HFB-27 EDFs, the drip points are -300 whereas again 
the highest value 4 t | np =306 there is for the SKM* EDF. 

The ground state deformation parameter fh of the Th 
chain and the corresponding deformation energy are dis¬ 
played in Figs. [5] and [6] correspondingly. Again all five 
EDFs lead to similar results until A =* 280, and there 
is rather extended region of spherical nuclei around this 
point. For all EDFs under consideration, the positive de¬ 
formations appear in the drip region, the values of defor¬ 
mation parameter reaching fh — 0.2. 
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Figure 5. /? 2 values in the thorium chain for different EDFs. Figure 8. /3 2 values in the plutonium chain for different EDFs. 



Figure 6. £ def values in the thorium chain for different EDFs. 



Figure 7. S 2 „ values in the plutonium chain for different EDFs. 


Lastly, we investigate the plutonium isotopic chain. 
The two-neutron separation energies are shown in Fig. |7] 
Qualitatively, the picture reminds very much that for the 
U case, see Fig. Q] Again, all EDFs except SKM* repro¬ 
duce reasonably the experimental data, and for the SKM*, 
the two-neutron separation energies are higher, especially 
in the drip region. Again the corresponding drip line point 
value ,4 ( i np (SKM*)=324 is much higher that those for other 

EDljShe deformation parameter /L and the corresponding 
deformation energy £def(/? 2 ) are shown in Figs. [8] and [9] 
correspondingly. The HFB-17 and HFB-27 EDF reveal a 
region of oblate deformations at 266 < A < 268. None of 
other EDFs confirm it. Again, this disagreement is a con¬ 


sequence of a competition of two close by energy minima, 
as was discussed earlier. For HFB-17 and HFB-27 EDFs, 
the oblate minimum has lower energy, whereas for other 
EDFs, the prolate minimum has a lower energy. Similarly 
as in the uranium chain, the spherical region ends for the 
FaNDF 0 EDF a bit earlier than for others. One more pe¬ 
culiarity occurs for the Fayans EDF: after the last spheri¬ 
cal nucleus 28() Pu, a small oblate deformation fh- — 0.05 
appears in 2X2 Pu, which than becomes prolate deformed, 
/?2-0.1 in 284 Pu. In this case, there is a competition of three 
close by minima, a spherical one and two deformed ones. 
The spherical minimum is lowest in 28() Pu, the oblate one 
in 2S2 Pu and the prolate minimum in 284 Pu. Note that such 
a small negative deformation in 282 Pu is almost invisible 
in the deformation energy curve, in Fig. [9] Similarly as in 
the uranium chain, all the EDFs under consideration lead 
to a strong prolate deformations in the drip line region, and 
again the FaNDF 0 deformation is the largest one. 


3 Two-neutron drip line 


In this section, we analyze the two-neutron drip line for a 
set of elements with an even Z value, from Pb up to Fm. 
These nuclei are important for analysis of the r-process in 
stars [35-37]. As in the previous Section, we compare pre¬ 
dictions from the FaNDF 0 functional with those from four 
Skyrme EDFs. The results are presented in Table 1. The 
values in the SLy4 and SKM* columns are found using 
the code [132]. They coincide or are very close to the cor¬ 
responding results presented in [33]. Small differences can 
be explained with some distinctions in the calculation de¬ 
tails, for example by the number of used oscillator shells. 
The values in HFB-17 and HFB-27 are taken from fToh . 
For each A," p , the corresponding value of the deforma¬ 
tion parameter is given in brackets. We see that the major 
part of the EDFs under consideration predict a spherical 
form for the drip line nuclei in the region from Pb until Ra, 
whereas all of them, from Th and heavier, are deformed. 

As a rule, predictions from the SKM* EDF are sig¬ 
nificantly higher compared to others. The difference is 
especially large for Rn and Ra elements for which this 
EDF predicts a strong deformation in the drip line re¬ 
gion, whereas these nuclei remain to be spherical for other 
EDFs. For the SKM* EDF, a strong irregularity in the 
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Table 1. Two-neutron drip line points A* lp (Z), from Pb to Fm, for different EDFs. In addition, the deformation /?2 of the drip line 

nucleus is given in brackets. 


element 

FaNDF 0 

SLy4 

SkM* 

HFB-17 

HFB-27 

Pb 

266 [0.00] 

266 [0.00] 

272 [0.00] 

266 [0.00] 

266 [0.00] 

Po 

270 [0.00] 

272 [0.00] 

280 [0.00] 

268 [0.00] 

268 [0.00] 

Rn 

268 [0.00] 

274 [0.00] 

298 [0.20] 

272 [0.00] 

274 [0.00] 

Ra 

272 [0.00] 

276 [0.00] 

304 [0.21] 

>286 

>286 

Th 

294 [0.24] 

292 [0.19] 

306 [0.21] 

300 [0.23] 

298 [0.22] 

U 

300 [0.26] 

298 [0.22] 

318 [0.21] 

308 [0.25] 

304 [0.24] 

Pu 

308 [0.27] 

304 [0.24] 

324 [0.20] 

314 [0.25] 

310 [0.25] 

Cm 

312 [0.28] 

310 [0.25] 

336 [0.14] 

320 [0.24] 

316 [0.24] 

Cf 

318 [0.27] 

316 [0.25] 

354 [0.00] 

324 [0.24] 

322 [0.24] 

Fm 

324 [0.26] 

322 [0.24] 

342 [0.23] 

330 [0.23] 

328 [0.23] 



Figure 9. £def values in the plutonium chain for different EDFs. 


^drf p (Z) dependence occurs for Cf. The anomalously high 
value of A^‘ P (SKM*) is explained with the the competi¬ 
tion between the spherical and prolate minima of the total 
energy, the spherical one being lower, whereas the pro¬ 
late ones are lower in all the neighboring elements. The 
predictions from the HFB-17 and HFB-27 EDFs are, as a 
rule, rather close to each other. The maximal difference 
between the corresponding values of A*‘ p equals to 4. As 
to the FaNDF 0 predictions, they turned out to be very close 
to those of the SFy4 EDF. The difference between them 
and both HFB EDFs is also quite moderate. 


A remark should be made concerning the relation be¬ 
tween A* lp values and one-neutron drip line points A^ np . 
Usually, due to the pairing effect, the inequality Ajj np <A A ^ 
is valid 123 El- 


Finally, we would like to stress, that in particularly 
such a heavy region of the nuclear chart, extrapolation of 
current EDF models up to the drip line is prone to large 
uncertainties. Therefore, the results presented here should 
be taken with a typical uncertainty bar of +10 mass units, 
or more 113 311 . An additional shift of the two-neutron drip 
line may occur due to the effect of p-dcpendcnce of the 
EDF parameters Hii, as discussed in the Introduction. 


4 Conclusions 


Fayans energy density functional FaNDF 0 is applied to 
nuclei around the uranium region. For Th, U, and Pu 
isotopic chains, the two-neutron separation energies S’ 
the ground state quadrupole deformation parameter /T 
and the corresponding deformation energies F’dcK/T) are 
found and compared with predictions from several Skyrme 
EDFs. Those from the SLy4 and SkM* EDFs were cal¬ 
culated by using the code fl32ll . whereas the HFB-17 and 
HFB-27 predictions are taken from Ei. For the major 
part of nuclei with known experimental S 2 n values, the re¬ 
sults obtained with the Fayans and SLy4 EDFs are rather 
close to those with two HFB EDFs, the HFB-17 one be¬ 
ing highly accurate in the overall description of nuclear 
masses |9|]. The SKM* EDF overestimates S 2 ,, values 
leading for these three elements to the two-neutron drip 
line point A," p values noticeably higher than those ob¬ 
tained with the other EDFs under consideration. The de¬ 
formation characteristics, the deformation parameter /B 
and the corresponding deformation energy £def(/? 2 X for the 
FaNDF 0 EDF, are also in the overall agreement with those 
of SLy4 and two the HFB EDFs. The SKM* predictions 
are again different, especially around the drip line region. 


The two-neutron drip line is found for FaNDF 0 , SLy4 
and SkM* EDFs for a set of elements with even proton 
number, from Pb up to Fm. This part of the nuclear chart 


is important for the study of the r-process in stars Q35 


37]. The consideration is made within a standard approach 


with fixed EDF parameters found mainly for stable nuclei. 
Within such approach, there is an overall agreement be¬ 
tween all EDFs under consideration except SKM*, the lat¬ 
ter predicting sufficiently higher A^ llp values. 


To conclude, the ground state properties of deformed 
nuclei in the uranium region, predicted from the FaNDF 0 
EDF, are found to be rather similar to those from several 
popular Skyrme EDFs. The same also holds for the predic¬ 
tion of the two-neutron drip line for even proton number 
elements from Pb to Fm. With such estimates, however, 
one should be careful since they contain a lot of uncertain¬ 
ties. Nevertheless, the FaNDF 0 prediction for the drip line 
seems to be in the line with typical Skyrme EDFs. 
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